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The pseudo-rapidity asymmetry and centrality dependence of charged hadron spectra in d+Au 
collisions at ^/snn = 200 GeV are presented. The charged particle density at mid-rapidity, its 
pseudorapidity asymmetry and centrality dependence are reasonably reproduced by a Multi-Phase 
Transport model, by HIJING, and by the latest calculations in a saturation model. Ratios of 
transverse momentum spectra between backward and forward pseudorapidity are above unity for 
Pt below 5 GeV/c. The ratio of central to peripheral spectra in d+Au collisions shows enhancement 
at 2 < pr< 6 GeV/c, with a larger effect at backward rapidity than forward rapidity. Our measure- 
ments are in qualitative agreement with gluon saturation and in contrast to calculations based on 
incoherent multiple partonic scatterings. 
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Soft and hard scattering processes have distinctive ra- 
pidity and centrality dependences in the context of par- 
ticle production in d(p)+Au collisions. Models based on 
the Color Glass Condensate [HQ, HIJING 0, and Multi- 
Phase Transport (AMPT) |4| predict specific pseudora- 
pidity (77) and centrality dependence of produced particle 
density which can be directly compared to experimental 
measurements. The Cronin effect |5J — the enhancement 
of particle yield at intermediate transverse momentum 
(pt) with respect to binary collision scaling— has also 
been observed in d+Au collisions at RHIC [111 II II Gil 
For partonic processes such as the dominant g + g and 
q + g scatterings, the particle rapidity distribution can be 
evaluated in a pQCD-inspired framework that depends 
on the parton distribution functions and the underlying 
dynamics. For example, calculations of the Cronin effect 
based on incoherent initial multiple partonic scatterings 
and independent fragmentation |3j predict a unique ra- 
pidity asymmetry of particle production in d+Au colli- 
sions, where the backward-to-forward (negative rapidity 
(Au) to positive rapidity (d)) particle ratio is greater than 
unity at low px, goes below unity at intermediate pr, and 
approaches unity again at high px- The amplitude of the 
theoretical backward-to-forward particle ratios depends 
on the nuclear shadowing U ■ Calculations of shado wing 
alone, based on Regge theory and hard diffraction \U\ . 
are fairly successful in describing the observed suppres- 
sion of particle production at forward rapidity in d+Au 
collisions |l2j. The calculation in Ref. |l2j considers the 
spatial dependence of the shadowing, leading to an im- 
pact parameter dependence that goes beyond the simple 
geometrical scaling. Calculations in a gluon saturation 
model 0| predict a backward-to-forward particle ratio 
that is opposite to the predictions based on incoherent 
multiple partonic scatterings. In this approach, the par- 
ticle production is related to the high gluon density in 
the nucleus (nucleon). The asymmetry is greater than 
unity in the range of transverse momenta determined by 
the values of the saturation scale Q s (y) and the geomet- 
rical scale Q 2 s {y)/Q s ,min, where Q s ,min is at the onset of 
the gluon saturation. Recently, the quark recombination 
model was used to explain the Cronin effect as a final 
state effect 0] , implying a backward-to-forward particle 
ratio markedly different from that of the QCD-inspired 
formulation in !3j a nd similar to the predictions by a sat- 
uration model 13] . In this approach, the enhancement 
of particle production at intermediate pt is an extension 
from low pt due to the thermal parton and shower parton 
recombination 14] . 

The suppression of high transverse momentum par- 
ticles in central Au+Au collisions at RHIC can be de- 
scribed by both final state and initial state effects, such 
as jet quenching calculations that assume parton energy 



loss via gluon bremsstrahlung 0, 0] or gluon satura- 
tion 0. The measurement of particle production at 
mid-rapidity from d+Au collisions at RHIC 0, US 
favors the scenario that the suppression of high-pT par- 
ticles is primarily due to the final state interactions, i.e., 
processes after the hard partonic scattering. The quanti- 
tative features of high-p^ particle production in Au+Au 
collisions can be described by models that incorporate a 
combination of physical effects such as the Cronin effect, 
nuclear shadowing 0] , and parton energy loss 0, 0] . 
The Cronin effect and shadowing can be investigated 
in d(p)+Au collisions. The magnitude of these nuclear 
effects on particle production has a geometrical depen- 
dence due to the nuclear density distribution. The parti- 
cle production in d(p)+Au collisions at different rapidi- 
ties also reflects the dynamics of nuclear and Bjorken-x 
dependence of these effects. Therefore, the centrality, 
pseudorapidity and pr dependence of particle produc- 
tion in d(p)+Au collisions provides an essential baseline 
for understanding the underlying phenomena in Au+Au 
collisions. 

We present inclusive px spectra of charged hadrons 
over an 77 range of — l(Au-side) to +l(d-side) in d+Au 
collisions at ^/snn = 200 GeV with several collision cen- 
trality selections. For these measurements, the STAR 
Time-Projection Chamber (TPC) provided tracking 
of charged hadrons. The minimum bias trigger was de- 
fined by requiring that at least one beam-rapidity neu- 
tron impinge on the Zero Degree Calorimeter |2(1 | in the 
Au beam direction. The measured minimum bias cross 
section amounts to 95+3% of the total d+Au geometric 
cross section. Charged particle multiplicity within —3.8 
< 7] < —2.8 was measured by the Forward TPC [2l| 
in the Au beam direction and served as the basis for 
our d+Au centrality tagging scheme, as described in [|J. 
The d+Au centrality definition consists of three event 
centrality classes; the 0-20, 20-40 and 40-100 percentiles 
of the total d+Au cross section. A separate centrality 
tag, which requires that a single neutron impinge on the 
Zero Degree Calorimeter in the deuteron beam direction 
(ZDC-d), was also used. Our analysis was restricted to 
events with a primary vertex within 50 cm of the center 
of the TPC along the beam direction. This yielded a data 
set of 9.5 x 10 6 minimum bias events. Only tracks (with 
at least 15 measured points) with a projected distance of 
closest approach to the event primary vertex of less than 
3 cm were used in the analysis. 

Acceptance and TPC tracking efficiency corrections in 
various pseudorapidity regions and centrality classes were 
obtained by embedding simulated data into a real data 
sample. In the region of |?7| < 0.5, the tracking efficiency 
and acceptance above pr = 2.0 GeV/c were observed to 
reach a plateau of about 90% for all centrality classes. Ef- 
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FIG. 1: The pr spectra of charged hadrons. From the top, 
the open circles correspond to the 0-20%, 20-40%, minimum 
bias, and 40-100% centralities in —1.0 < r\ < —0.5. Similarly, 
the solid triangles, open squares, and solid squares correspond 
to pr spectra in 0.0 < Tj < 0.5, —0.5 < rj < 0.0, and 0.5 < 77 
< 1.0, respectively. Spectra have been scaled by the factors 
indicated in the figure. 



ficiency corrections using filtered HIJING |22J — HIJING 
events in a GEANT simulation of the detector — were also 
used; a maximum difference between HIJING and em- 
bedded data of about 3% was observed. Background due 
to weak decay products was accounted for using filtered 
HIJING. For the 0-20% most central events, the contam- 
inating signals are estimated at less than 18% for pt<1-0 
GeV/c, and for the 40-100% most peripheral events this 
was observed to be less than 12%. The background ex- 
ponentially decreases, and above pr = 1.0 GeV/c, the 
background is approximately 4%, exhibiting no strong 
dependence on centrality or pseudorapidity. A net uncer- 
tainty of 6% in the analysis corrections was determined 
by adding the efficiency and background correction un- 
certainties in quadrature. 

The transverse momentum spectra of primary charged 
hadrons for various pseudorapidity regions are shown in 
Fig. □ for the 0-20%, 20-40%, 40-100% centrality selec- 
tions, and for minimum bias events. In the region of 0.2 
< pr< 2.0 GeV/c, the charged hadron spectra were fitted 



30 
25 



20 u^...,. 



0-20% 

20-40% 

40-100% 

Minbias 

AMPT 

Saturation 




FIG. 2: (Color online) The pseudorapidity dependence of 
charged particle densities for various centrality classes. Par- 
ticle tracking efficiency and background corrections were car- 
ried out for each pseudorapidity bin (Ar/=0.1). The point-to- 
point systematic uncertainties shown for each distribution (in- 
dicated by bands) are the quadratic sum of the efficiency and 
background correction uncertainties; statistical uncertainties 
are negligible. The results of AMPT (with default parame- 
ters) and parton saturation are indicated by the dashed and 
solid lines, respectively. 



with a power law function: 
d 2 N 



prdprdr] (1+pr/po)™ 



(1) 



The integrated charged hadron multiplicity per unit of 
pseudorapidity dN/dr] was obtained by summing up the 
measured yields in the covered momentum range and us- 
ing the power-law function for extrapolation to px = 
GeV/c. Fig. |21 shows the pseudorapidity dependence of 
charged particle densities for various centrality classes. 
Calculations based on the ideas of gluon saturation Q 
in the Color Glass Condensate as well as the predictions 
of AMPT Q are also shown. Both models predict a 
similar pseudorapidity dependence of particle yields. It 
should be noted that the pseudorapidity and centrality 
dependence of charged particle yields generated by HI- 
JING |22J (without shadowing) are nearly identical to 
the AMPT results at mid-rapidity. There is a clear in- 
crease in the asymmetry of charged particle densities as a 
function of increasing centrality: a prominent pseudora- 
pidity dependence is observed for the 0-20% most central 
collisions, while peripheral collisions between gold nuclei 
and deuterons are akin to symmetric p+p collisions. The 
predictions of the gluon saturation model and AMPT are 
in good overall agreement with the data. 

We define a measured asymmetry by taking ratios of 
inclusive backward (Au-side) to forward (d-side) px spec- 
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FIG. 3: (Color online) The ratio of charged hadron spectra 
in the backward rapidity to forward rapidity region for min- 
imum bias and ZDC-ri neutron-tagged events. Calculations 
based on pQCD 3] (y=—l/y=l) for minimum bias events are 
also shown for cases with no shadowing (solid curve), HIJING 
shadowing (dashed curve), and EKS shadowing (dot-dashed 
curve) . Calculations in a gluon saturation model [l3l| for min- 
imum bias events are shown for 0.5 < \rj\ < 1.0 (filled circles 
with solid line) and for 0.0 < |ry| < 0.5 (open squares with 
solid line). 
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tra. Fig. [3] shows the pr dependence of the asymme- 
try for minimum bias and ZDC-ri neutron-tagged events. 
The ratio was taken between the — 1.0< rj <— 0.5 and 
0.5< 7} <1.0 as well as -0.5< r\ <0.0 and 0.0< r\ <0.5 
regions. An overall systematic uncertainty (indicated by 
the band) of less than 3% was assessed by taking the 
corresponding ratios between inclusive spectra measured 
by STAR in p+p collisions at the same energy, where 
an asymmetry is not expected to be present. The ratio 
taken within |ry| < 0.5 is nearly constant in px, with a 
maximum value of approximately 1.075. This indicates 
that there is a small disparity between the forward and 
backward regions immediately around r\ = 0. The ra- 
tio taken at higher pseudorapidity slowly increases with 
Pt up to about pr = 2.5 GeV/c, attaining a value of 
approximately 1.25. The ratio taken at higher pseudora- 
pidity approaches unity beyond pr = 5 GeV/c, indicating 
the absence of nuclear effects at high pt- For the ZDC- 
ri neutron-tagged events, the ratio exhibits nearly the 
same pt dependence as minimum bias events. Fig. 0^, 
illustrates the centrality dependence of the asymmetry 
in the region of 0.5 < \rj\ < 1.0. The asymmetry be- 
comes more prominent with increasing centrality, reach- 
ing a factor of about 1.35 for the most central events. 
The asymmetry in the region of 0.0 < \rj\ < 0.5, shown in 
Fig. 3t>, does not exhibit a strong centrality and px de- 
pendence. The neutron-tagged events have an average 
number of binary collisions, {Nun) = 2.9±0.2, well below 
the {Nbi n ) = 7.5 ±0.4 of the minimum bias data set. The 
events in which a single nucleon from the deuteron inter- 
acted with the Au nucleus comprise approximately half 



FIG. 4: (Color online) a) The centrality dependence of the ra- 
tio of charged hadron spectra in backward rapidity to forward 
rapidity (0.5 < |ry| < 1.0). The gluon saturation model cal- 
culations are also shown for the 0-20% (solid curve), 20-40% 
(dashed curve), and 40-100% (dot-dashed curve) centrality 
classes, b) The centrality dependence of the ratio of charged 
hadron spectra in backward rapidity to forward rapidity (0.0 
< ry| < 0.5). The gluon saturation model calculations are also 
shown for the 0-20% (solid curve), 20-40% (dashed curve), and 
40-100% (dot-dashed curve) centrality classes. 



of the 40-100% peripheral centrality class However, 
Fig. 13 shows that the 77 asymmetry ratios for minimum 
bias and neutron-tagged events are nearly identical. 

Particle production at mid-rapidity in d+Au collisions 
may include contributions from deuteron-side partons 
that have experienced multiple scatterings while travers- 
ing the gold nucleus, and from gold-side partons that may 
have been modified by nuclear effects. Also shown in 
Figgis the calculation of the asymmetry in the incoher- 
ent multiple partonic scattering framework with various 
nuclear shadowing parametcrizations: no nuclear shad- 
owing, the HIJING shadowing [23j], and the EKS shad- 
owing 24] parameterizations. The ratio, taken for mini- 
mum bias spectra at y = — 1 and y = 1, is below unity at 
Pt ~ 3-4 GeV/c and is a consequence of the increase in 
Pt for partons from the deuteron hemisphere. Our mea- 
surements disagree with the theoretical calculations Q 
and thus suggest that incoherent multiple scattering of 
partons in the initial state alone cannot reproduce the 
observed pseudorapidity asymmetry in the intermediate 
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Pt region. By the same token, the class of models that 
incorporate initial parton scattering [IJ 0, |25|, though 
capable of reproducing integrated observables such as 
charged particle yield asymmetries, may not adequately 
reproduce the px dependence of the asymmetry. In this 
respect, the pt dependence of the pseudorapidity asym- 
metry as illustrated by the backward-to-forward ratio of 
charged hadron spectra can serve as an important dis- 
criminator between models. 

The minimum bias gluon saturation results for the 
backward-to-forward ratio of charged hadron spectra, 
also shown in Fig. were obtained by performing a cal- 
culation identical to the one in Ref. |l3| on the basis of 
the method developed in 0, |2(| . In this approach, the 
asymmetry is greater than unity in the range of trans- 
verse momenta determined by the values of the satura- 
tion scale Q s (y) and the geometrical scale Ql{y) / Q s ,min- 
The calculated particle yield aysmmetry, evaluated over 
the same pseudorapidity range as the data, is in qual- 
itative agreement with our observations. The theoreti- 
cal asymmetry exhibits a stronger pt dependence than 
actually observed, overpredicting the magnitude of the 
asymmetry at high pseudorapidities. The centrality de- 
pendence of the backward-to-forward particle yields in 
a saturation model, illustrated in Fig. |2Ji and Fig. f4J), 
qualitatively reproduces the observed centrality depen- 
dence. Although the model calculations fail to describe 
the data in detail, they show the same trend of increas- 
ing asymmetry with increasin g ce ntrality. We note that 
some conventional models [lM l27j are able to reproduce 
the suppression of particle production at forward rapid- 
ity in d+Au collisions, which was thought to be a unique 
feature of gluon saturation 0, El Hij . It will be inter- 
esting to quantitatively compare our measurements with 
those calculations in the future. 

It should be noted that a strong particle dependence 
in the nuclear modification factor has been observed in 
this intermediat e pt region in both Au+Au [2£| and 
d+Au collisions [lfj| . Collective partonic effects at the 
hadron formation epoch such as parton coalescence or 
recombination yfj, bill |32, |33| have been proposed to ex- 
plain Au+Au results. The pseudorapidity asymmetry 
approaches unity at a px scale above 5 GeV/c, approx- 
imately the same pt scale above which the particle de- 
pendence of the nuclear modification factor disappears. 
The idea of recombination was modified to explain the 
Cronin effect and its particle dependence Q as a final 
state effect. In this approach, the enhancement of par- 
ticle production at intermediate pt is an extension from 
low pt due to the thermal parton and shower parton re- 
combination [T^|. qualitatively consistent with the mea- 
surements of the pseudorapidity asymmetry as a func- 
tion of pt- We should emphasize that the pseudorapidity 
asymmetry is not likely to be solely due to the change of 
particle composition. In the recombination model, the 
shower and thermal parton recombination not only en- 
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FIG. 5: (Color online) The ratio of central (0-20%) to periph- 
eral (40-100%) spectra in d+Au collisions for various pseudo- 
rapidity regions and in Au+Au collisions at mid-rapidity. 



hances the baryon production, but also the meson pro- 
duction ^4| . The pseudorapidity asymmetry of identified 
pion spectra and its quantitative comparison to models 
are important for further understanding of particle pro- 
duction at intermediate pt- 

Of similar interest is the ratio of d+Au central to pe- 
ripheral inclusive spectra 



pdAu _ 
n CP — 



(d 2 N/dp T dr,/(N bin ))\ 



central 



(d 2 N/dp T d V /(N hia ))\ 



pcriph 



(2) 



where d? N / dprdrj is the differential yield per event in 
collisions for a given centrality class and (N^i n ) is the 
mean numbers of binary collisions corresponding to this 
centrality. Using a Monte Carlo Glauber calculation, as 
described in [||, the mean number of binary collisions 
for the 0-20% and 40-100% centrality classes were deter- 
mined to be 15.0+1.1 and 4.0+0.3, respectively. Fig. 
shows the ratio of the central to peripheral spectra in 
d+Au collisions for various pseudorapidity regions. The 
error bars on each distribution are the quadratic sum 
of statistical and systematic uncertainties; the latter are 
due to uncertainties in our background subtraction tech- 
nique. An overall error of about 10% due to the uncer- 
tainty in normalization is indicated by the band on the 
left portion of the figure. The Rcp in Au+Au collisions 
at y/s^N = 200 GeV is shown on the bottom of the 
plot. i?pp u distributions for each pseudorapidity selec- 
tion exhibit a rise with increasing pt, exceeding unity 
at p T ~l-2 GeV/c. At low p T , the i?^p u distribution is 
highest for the most backward pseudorapidity region and 
systematically decreases the more forward in pseudora- 
pidity the ratio is taken. The trend in the pseudorapid- 
ity dependence indicates that the Cronin effect is more 
pronounced in the gold hemisphere of the collision, con- 
sistent with the measured asymmetry between backward 
and forward rapidity. Our measurement of i?pp u shows 
no significant suppression at pt of 2-6 GeV/c. This result 
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stands in contrast to the Au+Au measurements, where 
Rcp was observed to be well below unity for pt < 12 
GeV/c. The results for R$fp are consistent with calcu- 
lations in pQCD models incorporating both Cronin en- 
hancement and nuclear shadowing [23, HE |3|| 113, . 
However, the models based on incoherent parton scatter- 
ing at the initial stage fail to reproduce the rapidity de- 
pendence in both backward-to- forward ratios and R^p" ■ 

In summary, we have studied the centrality and pseu- 
dorapidity dependence of charged hadron production in 
d+Au collisions at ^/snn — 200 GeV. The inclusive 
charged hadron multiplicity is observed to be higher 
in the gold hemisphere than the deuteron hemisphere 
of the collision. The gluon saturation, HIJING, and 
AMPT models cannot be ruled out from the integrated 
charged particle pseudorapidity distributions. Ratios of 
backward-to-forward pseudorapidity transverse momen- 
tum distributions are above unity for px below 5 GeV/c. 
Our measurement of R(fp shows no suppression at pt of 
2-6 GeV/c, with the ratio taken at backward pseudora- 
pidities being slightly higher than at forward pseudora- 
pidities. The incoherent multiple scattering of partons 
in the initial state alone cannot reproduce the observed 
pseudorapidity asymmetry, while the latest calculations 
in a gluon saturation model stand in qualitative agree- 
ment with our observations. 
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